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Abstract 

Quantum molecular dynamic simulations are introduced to study the dynamical, electrical, and 
optical properties of carbon dioxide under dynamic compressions. The principal Hugoniot derived 
from the calculated equation of states is demonstrated to be well accordant with experimental re- 
sults. Molecular dissociation and recombination are investigated through pair correlation functions, 
and decomposition of carbon dioxide is found to be between 40 and 50 GPa along the Hugoniot, 
where nonmetal-metal transition is observed. In addition, the optical properties of shock com- 
pressed carbon dioxide are also theoretically predicted along the Hugoniot. 
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I. INTRODUCTION 

Knowledge of pressure-induced transformations for materials under extreme conditions, 
which requires accurate understandings of the thermophysical properties into new and com- 
plex region, has gained much interest in a large number of scientific and technological do- 
mains l|. Carbon dioxide (CO2) is one kind of the most important molecules in planetary 
science and chemical (or explosive) process. For instance, the magnetic fields and evolution 
of planets are closely related to the equation of states (EOS) and electrical properties of 
planetary interiors J2|, yj, which are composed of species containing C, H, N, O, etc.; Chem- 
ical reactions of carbonaceous molecules (such as C0 2 , CO), nitrogen (N 2 ), oxygen (0 2 ) 
and their mixture at high densities and temperatures are focused due to their presence in 
reacted chemical explosives |4H7|. F rom these aspects, clearly, an accurate description of 
thermodynamical properties of C0 2 at high pressures and temperatures is fundamental and 
proves itself to be indispensabe for developing realistic models of planetary interiors and 
reacted explosives. 

Solid phase diagram of carbon dioxide up to 80 GPa and molecular-nonmolecular tran- 
sition have already been probed through diamond-anvil cells and confocal micro-Raman 
spectroscopy [a-lll]. Solid specimens were reported to be stable until the formation of a- 



carbonia (an amorphous extended phase), which can only be revealed around 50 GPa at 
elevated temperature or above 65 GPa at ambient temperature [11]. Here, we concentrate 
on liquid specimens, which were first single shocked to 29 GPa and double-shocked to 56 



GPa by Schott [12J, then single shocked to 71 GPa by Nellis et al. 



131 ] using a two-stage 



light-gas gun. Theoretically, intermolecular potential method was applied to study the EOS 
of C0 2 under dynamic compressions, where the onset of decomposition was reported to be 
around 30 GPa (4500 K, 17 cm 3 /mol) along the Hugoniot [13J . Despite that the predicted 
EOS were in accord with experimental data, the electronic structure, which has been proved 
to be important in determining the dynamical, electrical, and optical properties of molecular 



14 



16], is still lacking due to the intrinsic approximations 



fluids under extreme conditions 
of this method. 

On the other hand, quantum molecular dynamic (QMD) simulations, where quantum- 
mechanical treatments are executed by combining classical molecular dynamics for the ions 
and density functional theory (DFT) for electrons, have already been proved to be successful 



in describing the thermodynamical properties of mono atomic molecules (He) 17j|, diatomic 



molecules (H 2 , N 2 , O2) |18|-|20|. and hydrocarbons (CeH 6 ) [l6| under extreme conditions. 
A fully quantum-mechanical description of CO2 under shock compressions is highly recom- 
mended to be presented and understood because of the following aspects: (i) C0 2 undergos 
dramatic transformations under shock compressions, such as dissociation and recombination 
of molecules, during which carbonaceous molecules combined with oxygen (O2) might be 
produced, and electron spin polarization should be seriously considered at this stage due to 
the fact that spin-triplet state governs the ground state of oxygen [21] ; (ii) Electrical conduc- 
tivity of molecular fluid, where electronic structure is of dominance, can be greatly influenced 
by thermo-activated electronic states under extreme conditions; (iii) Temperature, which is 
an important parameter in experimental determination of EOS, is usually difficult to be 
measured because of the uncertainty in determining the optical intensity for ultraviolet part 
of the spectrum, and QMD simulations are powerful tools to provide predictions. 

In the present work, DFT-based QMD simulations have been used to investigate ther- 
mophysical properties of CO2 under extreme conditions. The EOS and pair correlation 
functions (PCF) are determined through QMD simulations at equilibrium. Dynamic con- 
ductivity ct(oj) is calculated by Kubo-Greenwood formula, from which the dc conductivity 
(adc) is determined. Then, the dielectric function e{uS) and reflectivity are extracted. The 
rest of this paper is organized as follows. The simulation details are briefly described in 
Sec. [Til The PCF, which is used to study the dissociation of CO2, and the Hugoniot curve 
are given in Sec. IIIIAt In Sec. IIIIBl nonmetal-metal transition and optical properties are 
discussed. Finally, we close our paper with a summary of our main results in Sec. [TV] 

II. COMPUTATIONAL METHOD 



The Vienna Ab-initio Simulation Package (VASP) [22l . l23l | , which was developed at the 
Technical University of Vienna, has been employed to perform simulations for carbon diox- 
ide. The elements of our calculations consist of a series of volume-fixed supercells includ- 
ing N atoms, which are repeated periodically throughout the space. By involving Born- 
Oppenheimer approximation, electrons are quantum mechanically treated through plane- 



wave, finite-temperature (FT) DFT 24|, where the electronic states are populated accord- 



ing to Fermi-Dirac distributions at temperature T e . The exchange-correlation functional 



is determined by generalized gradient approximation (GGA) with the parametrization of 



Perdew-Wang 91 25]. The ion-electron interactions are represented by a projector aug- 



mented wave (PAW) pseudopotential [26|. Isokinetic ensemble (NVT) is adopted in present 
simulations, where the ionic temperature Tj is controlled by Nose thermostat [27| . and the 
system is kept in local thermodynamical equilibrium by setting the electron (T e ) and ion 
(Tj) temperatures to be equal. Electron spin polarization has also been taken into account, 
due to the possible existence of oxygen molecules in shocked CO2. 

The plane-wave cutoff energy is selected to be 600 eV so that the pressure is converged 
within 5% accuracy. T point and 4x4x4 Monkhorst-Pack scheme k points are used to sample 
the Brillouin zone in molecular dynamics simulation and electronic structure calculation, 
respectively, because EOS (conductivity) can only be modified within 5% (15%) for the 
selection of higher number of k points. 81 atoms (27 CO2 molecules) are included in the 
cubic supercell. The densities selected in our simulations range from 13 to 37.55 cm 3 /mol 
and temperatures between 218 and 10000 K, which highlight the regime of the Hugoniot. 
All the dynamic simulations are lasted for 3 ~ 6 ps, and the time steps for the integrations 
of atomic motion are 0.5 ~ 1 f s according to different densities (temperatures). Then, the 
subsequent 1 ps simulations are used to calculate EOS as running averages. 

III. RESULTS AND DISCUSSION 

A. EQUATION OF STATE 

Accurate understandings of the electrical and optical properties of CO2 depend on a 
precise description of dynamical properties, such as EOS. The EOS have been examined 
theoretically through the Rankine-Hugoniot equations, which follow from conservation of 
mass, momentum, and energy across the front of shock waves. The equations describe the 
locus of points in (E, P, V^-space satisfying the relation as follows: 

E 1 -E = ±(P 1 + P )(V-V ), (1) 

P 1 - P = p u s u p , (2) 

K=V (1-VO> ( 3 ) 
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FIG. 1: (Color online) The P-V (left panel) and T-P (right panel) Hugoniot curves of shocked 
carbon dioxide. Inset is the (u s , u p ) diagram. Here, filled red circles denote the present QMD 



results, while the dashed and dotted lines denote the intermolecular potential results 13|] with and 
without accounting for moleculardissociation, respectively. The experimental data reported in 
Ref. [12|] (open circles) and Ref. [13| (open squares) are also shown for comparison. 

where E, P, V denote internal energy, pressure, volume, and subscripts and 1 present the 
initial and shocked state, respectively. In Eqs. (j2j) and (J3J), u s is the velocity of the shock 
wave and u p corresponds to the mass velocity of the material behind the shock front. In the 
present work, the initial density of C0 2 is p — 1.17 g/cm 3 (Vo=37.55 cm 3 /mol) and the liquid 
specimen is controlled at a temperature of 218 K, where the internal energy is Eq=— 22.95 
eV/molecule. The initial pressure Pq can be treated approximately as zero compared to the 
high pressure of shocked states along the Hugoniot. The Hugoniot points are obtained as 
follows: (i) smooth functions are used to fit the internal energy and pressure in terms of 
temperature at sampled density; (ii) then, Hugoniot points are derived from Eq. ([[]). The 
results are summarized in Table [H 

Comparisons between our simulated Hugoniot curve and results from experimental mea- 
surement and intermolecular potential method are shown in Fig. [1] (left panel), where our 



TABLE I: Hugoniot pressure and temperature points derived from QMD simulations at a series of 
molar volumes. 



V (cm 3 /mol) P (GPa) T (K) 



22 


5.72 


459 


21 


6.91 


568 


20 


10.05 


873 


19 


14.95 


1374 


18 


19.53 


1909 


17 


27.39 


2821 


16 


36.24 


3873 


15 


41.62 


4175 


14 


49.89 


4968 


13 


73.82 


7955 



results show good agreement with experiment. The QMD simulated pressure and temper- 
ature along the Hugoniot curve for carbon dioxide shows a systematic behavior, except for 
the region with a small break in slope around 40 to 50 GPa, where structural transitions 
characterize the EOS of molecular fluid. Although the P-V curve by intermolecular po- 
tential method is consistent with experiments, the predicted temperature remains too high 
compared with our calculations [see Fig. [T] (right panel)]. 

Dynamical properties and structural transitions of CO2 under extreme conditions have 
been examined through PCF, which represents the possibility of finding a particle at a 
distance r from a reference atom. Figure [2] shows PCF and charge density distribution (at 
equilibrium) of four different molar volume points, and the structural transitions of CO2 
are clearly reflected by PCF. At low pressures [Fig. EJ^a)], ideal molecular fluid is suggested 
by PCF distribution of C-O, which exhibits a main peak around 1.16 A. With the increase 
of pressure along the Hugoniot, this main peak in C-0 PCF reduces in amplitude and get 
broadened, which is attributed to the thermo-induced vibrations of C-0 bonds around the 
equilibrium length. Dissociation of CO2, which leads to soften behavior of the Hugoniot 
[see Fig. [1] (left panel)], has been found at about 40 to 50 GPa. It is indicated that carbon 
atoms tend to form connections between each other [Fig. GO^c)], then, with further increase 
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FIG. 2: (Color online) Calculated pair correlation function for C-0 (blue line), C-C (black line), 
0-0 (red line) at four densities of CO2 along the Hugoniot. The atomic structure, where carbon 
and oxygen atoms are denoted by gray and red balls respectively, and the relative iso-surface of 
charge density (blue regimes) are also provided in the insets, (a) V = 19 cm 3 /mol, T = 1374 K; 
(b) V = 16 cm 3 /niol, T = 3873 K; (c) V = 15 cm 3 /rnol, T = 4175 K; (d) V = 13 cm 3 /niol, T = 
7955 K. 

of pressure, hybrid C-C bonds (peak around 1.45 A, which lies between typical sp 3 and sp A 
hybrid C-C bond length), diatomic oxygen (peak around 1.23 A) and mono-atomic oxygen 
are formed in the shocked system [Fig. [21(d)]. Furthermore, we have examined the charge 
density difference between spin-up and spin-down electrons, which is as small as negligible, 
and the results suggest that electron spin polarization of the shock produced mono-atomic 
and diatomic oxygen are suppressed at this stage. 

B. DYNAMIC CONDUCTIVITY 



Great controversies have been raised since the detection of nonmetal-metal transition 
of such diatomic molecules as hydrogen and oxygen 



28 



29( in isentropic compressions. 



The links between nonmetal-metal transition and dissociation of molecules under dynamic 
compression are of particular significance. Here, we examine the dynamic conductivity of 
CO2 according to Kubo-Greenwood formula: 



N N 



°M = 5§s5>( k ) £ ££[/(*, k) - /(e;,k)] 

k j=l 8=1 a=l 

x|(^,k|V a |^, k )| 2 5(e Jik - e i>k - fuj), 



(4) 



where Q is the volume of the supercell. The i and j summations range over N discrete bands 
included in the calculation. The a sum is over the three spatial directions. /(e i? k) describes 
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FIG. 3: Calculated dc conductivity of shocked carbon dioxide along the Hugoniot. 

the occupation of the ith band, with the corresponding energy e^k and the wavefunction 
^k at k. u>(k) is the k-point weighting factor. 

Along the Hugoniot, dynamic conductivity have been calculated as the average of thirty 
atomic snapshots at equilibrium. The dc conductivity (<7dc), which follows the static limit 
(o; — ?> 0) of <Ji(u), is then extracted and plotted in Fig. [3j At pressures below 20 GPa, <Jd c 
can be neglected, and the results suggest an insulating molecular fluid. Thermal activation 
of conductivity is raised from 20 GPa, then metallization (a^c > 1000 n _1 cm _1 ) of CO2, 
which is accompanied with dissociation of molecules, is reached at 40 to 50 GPa. 

Emissivity of spectrum, which is closely related to optical reflectivity, is important in 
determining temperature in experiment. Optical reflectivity can be derived from dynamic 
conductivity. The imaginary part 0*2(0;) can be obtained from the Kramers-Kronig relation: 



0*2(0;] 






0\\y)us 

v 2 — u 2 



dv. 



(5) 



where P is the principal value of the integral. Then dielectric functions can be derived from 
the two parts of the conductivity: 



4-7T 

ei(w) = 1 0*2(0;), 

O" 



(6) 



4-7T 



£2(0; 



U) 



-0*1(0;). 



(7) 
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FIG. 4: (Color online) Optical reflectivity of shocked carbon dioxide for wavelengths of 404 (black 
line), 808 (red line), and 1064 nm (blue line) along the Hugoniot. 

The real part n(u) and the imaginary part k(u) of the refraction index have the following 
relations with dielectric function: 



n(u) = \l~[\e(u)\+e 1 (u)], 



(8) 



k{u) 



-[|c( w )|-ei(w)]. 



(9) 



The index of refraction is useful for evaluating optical properties such as the reflectivity r(u) 
and absorption coefficient a(uj): 

[1 - n(u)] 2 + k(u) 2 



r(ui) 



(10) 



[l + n(u)] 2 + k{u) 2 ' 

We have examined the optical reflectance at three wavelengths (404, 808, and 1064 nm) 
spanning the visible spectrum, and the results are shown along the Hugoniot in Fig. HI 
which clearly shows the pressure-induced change in reflectivity (from 0.05 to 0.35 ~ 0.55). 
For the pressure between 20 and 40 GPa, the increase in reflectance is retarded, then sharp 
increase exists around 40 to 50 GPa due to the nonmetal-metal transition. The results could 
be inspected by future experiments. 



IV. CONCLUSION 

In summary, we have studied the thermodynamical properties of carbon dioxide under 
extreme conditions. The EOS obtained from QMD simulations are used to determine the 
Hugoniot curve, which shows good agreement with available experimental data. Our cal- 
culations indicate that decomposition of carbon dioxide begins at 40 to 50 GPa along the 
Hugoniot, where carbon atoms intend to form connections between each other, then di- 
atomic and monoatomic oxygen are produced. The nonmetal- metal transition, as well as 
the change in optical reflectance, are also found in the present work. 
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